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DtSSPLA 


m 

Introduction 


This  document  is  intended  to  give  the  user  a  basic 
understanding  of^  Pi  sol  ay  A nteorated^ Software  System  and  Plotting 
Language.  Ail  examples  were  taken  from  implementations  at  David 
^Taylor  Naval  Ship  R  &  D  Center  to  show  actual  scientific  uses. 

All  examples  were  run,  unless  specifically  stated,  on  the 
VAXcluster  using  version  9*2  of  DISSPLA.  This  version  will  soon 
be  available  on  the  CYBERs  and  aside  from  certain  machine 
dependencies,  all  examples  will  then  run  on  the  CYBERs.  In 
interactive  examples,  underlined  commands  are  typed  by  the  user.  <T— 


Convent i ons 


Throughout  this  document  the  following  argument  naming 


conventions 

hold: 

ARGUMENTS: 

X.  IX,  LX 

Refer  to  the  X-axis 

Y,  IY,  LY 

- 

Refer  to  the  Y-axis 

1,  J,  K,  M,  N 

- 

Integer  Values 

All  Others 

- 

Real  Values 

L 

_ RAY 

- 

String  of  text 

Array  of  values 

NOTATION: 

{...} 

- 

Within  a  program  listing 
are  comments 

(under  1 i 

ne) 

- 

In  examples  are  user  type 

P/s 

parameter  setting 

[""Accession  For 

NTIS  GRA&I 
DTIC  TAB 
Unannounced 
Justif lcatlo 

% 

□ 

n  _____ 

•Rv 

Distr 

A\*n  i 

r,  *,  .-•  + 

UdL 

Lfcutlon 

Labilit 

Avail 

Spec 

/ 

y  Codes 
ind/or 

Lai 
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Devices 


The  call  to  a  device  MUST  be  the  first  call  in  a  DISSPLA 
program.  The  following  devices  are  currently  active;  Other 
devices  may  be  activated  at  the  user's  request,  these  calls  are 
supported  by  version  9*2  which  is  currently  only  on  the 
VAXcluster  and  should  soon  be  on  the  CYBERs. 


Tektronix  4000  Series: 


CALL  TK4006  ( 

call  moio  { 

CALL  TK4012  ( 
CALL  moi3  ( 
CALL  TK4014  ( 
CALL  TKU015  ( 

call  moi6  { 
CALL  TKU025 
CALL  TK4027 
CALL  TK14051  ( 
CALL  TK4052  ( 
CALL  TK4054  ( 

CALL  TEKALL  ( 


icps  ) 
icps  ) 
icps  ) 
icps  ) 

icps,  idevc  ) 
icps,  idevc  ) 
icps,  idevc  ) 


icps  ) 
icps  ) 

icps,  idevc  ) 
i model,  icps, 


0,  idevc,  0  ) 


CALL  PTEKAL  {will  prompt  for  all  information 


icps  -  Baud  rate  (30,  120,  240,  480,  960) 


idevc  -  is  1  for  high  resolution  (4096) 
-  is  0  for  low  resolution  (1024) 


Tektronix  4100  Series: 

CALL  TK41  (  IMODEL  ) 

IMODEL  -  4105,  4107,  4109,  4112,  4113,  4114,  4115,  4116 
CALL  PTK41  {will  prompt  for  model  number  } 
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DEC  Terminals: 


CALL  REGIS  (  ioptin,  idevc  ) 


i opt  in  -  1  (  VT125  ) 

-  2  (  VK 100  ) 

-  3  (  VT2*»0,  VT2M  ) 

-  1»  (DEC  PRO  350  ) 


idevc  -  0  (Monochrome  Monitor) 

-  1  (Color  Monitor  ) 

CALL  PREGIS  {will  prompt  for  information} 


CALCOMP: 


CALL  CALCOMP  (0,  0,  10) 

The  file  CALC0MP0UT.DAT  is  created  for  processing  to 
tape  using  the  procedure  VSYS:CALCD2T  on  the  VAX. 

The  file  TAPE10  is  created  for  processing  to  tape  on 
either  of  the  CYBERs. 


META  FILE: 


CALL  COMPRS 

The  file  DISPLOT.DAT  is  created  for  processing  by  the 
Post-Processor  on  the  VAX.  (RUN  VSYS-.DISPOP) 

The  file  PLFILE  is  created  for  processing  by  the 
Post-Processor  on  either  of  the  CYBERs.  (DISPOST) 


1 1  legal  Uni t  Numbers 

Unit  Numbers  used  by  DISSPLA  which  CANNOT  be  used  in  your 


program: 

31.32,33 

Scratch  files 

90,91,93 

Mapping  and  Landblan 

95 

META  file  for  COMPRS 

96 

Font  files 

97 

Scratch  file 

10 

CALCOMP  file 

DISSPU 
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This  is  an  example  of  running  a  DISSPLA  program  on  the 
VAXcluster . 


# 


PROGRAM  EXAMPLE! 
DIMENSION  X  00)  ,  Y  ( 1 0) 


{This  example  is  for  the  VAX  } 


DATA 

Y  /l, 

2,  3.  4.  5. 

DATA 

X  /l. 

2.  3.  4,  5. 

CALL 

COMPRS 

CALL 

BGNPL 

(0) 

CALL 

TITLE 

('This  is  a 

• 

1  Y-Ax i s$ ' , 

CALL 

GRAF 

(0.,  1.,  10 

CALL 

CURVE 

(X,  Y,  10, 

CALL 

ENDPL 

(0) 

CALL 

DONEPL 

END 

fortran  myf i le 
d i s 1 i nk  myf i le 

100,  6.,  8.) 

>.  o.,  1.,  10.) 


Other  Libraries  (y  or  n)  :  n 
$  run  myf i le 


PLOTTING  COMMENCING 


NO.  OF  FIRST  PLOT  0 


END  OF  DISSPLA  9-2  —  894  VECTORS  IN  1  PLOTS. 

RUN  ON  1/17/86  USING  SERIAL  NUMBER  3105  AT  DTNSRDC  VAX 
PROPRIETARY  SOFTWARE  PRODUCT  OF  1SSC0,  SAN  DIEGO,  CA. 

403  VIRTUAL  STORAGE  REFERENCES;  5  READS;  0  WRITES. 

$  run  vsys:d i spop 
ENTER  DEVICE  TYPE . 

1-TEK  40xx,  2-CALCOMP  1051,  3-TEK  4lxx,  4-REGIS  (VT240) ,  5-PRINTER  PLOT 


{  Enter  your  device  } 


8:46 


ENTER  MODEL  NUMBER 
4105 

ENTER  POST-PROCESSOR  DIRECTIVES 
<CR> 

PLOT  FILE  GENERATED  BY  Kevin  Brady 

AT  DTNSRDC  VAX 
ON  JAN  17,  1986 


4 

'*»i 

SI 


S' 

>a< 

I 
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This  command  procedure  on  the  VAXc luster  will  do  the  entire 
process  with  no  user  interaction: 

$  on  error  then  exit 
$  again: 

$  if  pi  .eqs.  ""  then  inquire  pi  "Filename" 

$  I  Keep  prompting  for  filename 
$  if  pi  .eqs.  ""  then  goto  again 
$  fortran  1  pi 1 

$  !  libraries  linked  twice  because  of  nowrap  search 

$  link  'pi  1 ,vsys:di sintf/I ibrary.disspla/l ibrary,- 
d i s i ntf /I ibrary.disspl a/ library 

$  run  'pi 1 

$  run  vsys:dispop 

3 

1*105  {  This  line  must  be  changed  to  your  device  } 

$  exit 


.‘i*| 


I 

ft 


1 


& 


I 
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The  ssme  program  run  interactively  on  either  of  the  CYBERs  to 
create  a  CALCOMP  plot  : 

PROGRAM  EXAMPLE  1 
DIMENSION  X  (10) .  Y (10) 


DATA 

Y  /l,  2 

.  3,  4.  5.  6, 

7.  8,  9,  10/ 

DATA 

X  /l,  2 

.  3.  4.  5.  6. 

7.  8,  9.  10/ 

CALL 

COMPRS 

CALL 

BGNPL 

(0) 

CALL 

HEIGHT 

(0.2) 

CALL 

TITLE 

("THIS  IS  A  SIMPLE  GRAPH",  22,  "X-AXIS",  6 

T- 

"Y-AXIS",  6, 

6..  8.) 

CALL 

GRAF 

(0 . ,  1 . ,  10., 

0.,  1.,  10.) 

CALL 

CURVE 

(X,  Y,  10,  1) 

CALL 

ENDPL 

(0) 

{Double  quotes  are  } 

CALL 

DONEP 

{required  on  CYBERs} 

END 

COMMAND-CONNECT. OUTPUT 
COMMAND-ATTACH. INF ILE.MYPLOT. ID-CAKB 
AT  CY-  001  SN-SYSSET 
C OMM AN D -ATTACH , D1SSPLA 

PFN  IS  DISSPLA 
AT  CY-  006  SN-SYSSET1 
COMMAND-REQUEST. TAPE  10. *PF 
COMMAND- ATTACH. NSRDC 

PFN  IS  NSRDC 
AT  CY*  105  SN-SYSSET1 
COMMAND-LIBRARY. DISSPLA. NSRDC 
COMMAND -FTN5 . 1-INF  I LE , L-0 

56100  CM  STORAGE  USED. 

0.020  CP  SECONDS  COMPILATION  TIME. 

COMMAND-LGO 


PLOTTING  COMMENCING 


.  DISSPLA  VERSION  8.2 

NO.  OF  FIRST  PLOT  0 


END  OF  DISSPLA  8.2  --  897  VECTORS  GENERATED  IN  1  PLOT  FRAMES. 

-ISSCO-  4186  SORRENTO  VALLEY  BLVD..SAN  DIEGO  CALIF.  92121 

DISSPLA  IS  A  CONFIDENTIAL  PROPRIETARY  PRODUCT  OF  ISSCO  AND  ITS  USE 
IS  SUBJECT  TO  A  NONDISSEMINATION  AND  NONDISCLOSURE  AGREEMENT. 
NON-cATAL  LOADER  ERRORS  -  SEE  MAP 
ENO  PLOTIT 

61300  MAXIMUM  EXECUTION  FL. 

0.132  CP  SECONDS  EXECUTION  TIME. 
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COMMAND-RETURN, D I SSPLA.NSRDC 
COMMAN D -ATTACH ,  DIJsPOST 

PFN  IS  DISPOST 
AT  CY*  007  SN=SYSSET1 
COMMAND-LIBRARY. DISPOST 
COMMAND-POP936 

**********  CALCOMP  POST-PROCESSOR  ********** 
INPUT  DIRECTIVES  {  Input  space  then  <CR>  } 


PLOT  FILE  GENERATED  BY  CAKB0E2  AT  10.05.27  ON  01/17/86 
1  THE  NUMBER  OF  CALCOMP  BLOCKS  IN  THIS  FILE  IS  3 

1  PLOTS  HAVE  BEEN  PROCESSED 

. .  END  OF  POSTPROCESSOR  . 

END  POP936 

COMMAND-CATALOG, TAPE10,  I  D*CAKB 
NEWCYCLE  CATALOG 
RP  *  030  DAYS 

CT  ID*  CAKB  PFN-TAPE10 

CT  CY*  002  00000016  PRUS  $0000.04  /DAY 

CT  SN=  SYSSET 


COMMAND- 
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The  following  procedure,  stored  in  the  user's  permanent  file 
PLOTIT,  will  let  him  view  his  output  on  a  Tektronix  terminal: 

.PROC.PLOTIT. 

ATTACH ,  I NF  I LE .MYPLOT, I D*CAKB . 

ATTACH, 01SSPLA. 

REQUEST, PLFILE,*PF. 

ATTACH, NSRDC. 

LIBRARY, DISSPLA.NSRDC. 

FTN5, l-INFILE.L-O. 

LGO. 

CATALOG, PLF I LE, ID-CAKB. 

RETURN, DISSPLA.NSRDC. 

ATTACH, D I SPOST. 

L I BRARY , D I SPOST. 

TEK300. 

REVERT. 

To  execute  type: 

C OMMAN D -ATTACH , PLOTIT , ID-xxxx 
COMMAND-PLOT  IT 


NOTE:  Currently  on  the  CYBERs  there  are  two  versions  of  the 
Tektronix  post-processor,  the  one  attached  above  is  used  for 
color  plots  and  WILL  NOT  stop  between  plots  if  there  are  more 
than  one.  To  view  multiple  plots  with  a  pause  between  each  one, 
the  other  version  of  the  post  processor  must  be  used.  Instead  of 
the  "ATTACH , D I SPOST"  and  "L I BRARY, D I SPOST" ,  MSFETCH  the  file 
TEK300  ("MSFETCH, TEK300, UN-CSYS")  and  continue. 
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Level  Structure 


The  axis  system  cannot  be  drawn  until  the  sub-plot  area  and 
physical  origin  are  defined.  Similarly,  curves  cannot  be  drawn 
until  the  axis  system  is  established.  Thus,  DISSPLA  has  a  level 
structure  to  ensure  that  all  necessary  information  is  present 
during  the  construction  of  a  plot.  An  error  message  is  printed  if 
a  routine  is  called  at  the  wrong  level. 

The  DISSPLA  levels  are: 

0  -  before  device  initialization 

1  -  after  device  initialization 


2  -  after  page  border,  physical  or 

igin,  and 

subplot 

area  defined 

3  -  after  axis  system  is  i 

def  i  ned 

Subrout i ne 

Level 

Level 

at 

after 

cal  1 

cal  1 

CALL 

DEVICE 

0 

1 

CALL 

AREA2D  (XAX1S,  YAXIS) 

1 

2 

CALL 

XNAME  (LXNAME ,  YXNAME) 

2 

p/s 

CALL 

YNAME  (LYNAME,  IYNAME) 

2 

p/s 

CALL 

HEADIN  (LHEAO,  1  HEAD ,  HTMULT 

,  NUNES) 

2,3 

same 

CALL 

GRAF  (XORG ,  XSTP,  XMAX,  YORG 

,  YSTP,  YMAX)  2 

3 

CALL 

CURVE  (XARAY ,  YARRAY,  NPTS, 

1  MARK) 

3 

same 

CALL 

ENDPL  ( 1  PLOT) 

2,3 

1 

CALL  DONEPL 


1 


0 
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Subrout i ne  Level  Level 

at  after 

cal  1  cal  1 

To  change  page  si2e  from  8.5  x  11  inch  default: 

CALL  PAGE  (XPAGE,  YPAGE)  1  p/s 

To  position  the  physical  origin: 

CALL  PHYSOR  (XPHYS,  YPHYS)  1  p/s 

To  re-position  the  origin  relative  to  PHYSOR: 

CALL  OREL  (XOREL,  YOREL)  1  p/s 

To  end  the  plot  and  remain  on  same  page: 

CALL  ENDGR  (IPLOT)  2.3  1 

To  frame  the  subplot  area: 

CALL  FRAME  2.3  same 


To  understand  the  concept  of  setting  up  the  plot,  consider 
the  following  example.  The  page  will  be  set  to  2V  X  11",  and 
will  contain  two  subplots.  The  lower  left  corner  of  the  first 
plot  is  positioned  at  the  point  (3", 2")  by  PHYSOR.  AREA2D  then 
sets  the  current  plotting  area  to  8"  X  6".  Now  the  plotting  can 
begin  (although  no  plots  are  actually  done  in  this  example).  The 
first  plot  is  then  ended  by  ENDGR,  which  tells  DISSPLA  you  are 
done  with  the  first  plot  but  will  rema i n  on  the  same  page  to  do 
another  plot.  The  lower  left  corner  of  the  second  plot  is  then 
put  at  (13". 2")  by  OREL  (moved  relative  to  previous  origin). 
AREA2D  is  called  again  to  define  the  second  plot  area,  and 
plotting  could  commence.  This  plot  is  ended  by  ENDGR  again  then 
ENDPL  terminates  the  plot  on  this  page. 
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PROGRAM  EXAMPLE2 


CALL  COMPRS 


CALL  BGNPL  (0) 

CALL  PAGE  (24.,  11.) 
CALL  PHYSOR  (3-.  2.) 
CALL  AREA20  (8.,  6.) 


CALL  FRAME 


{  Level  0  to  Level  1  } 


{  Level  1  to  Level  2  } 


CALL  HEADIN  ('THIS  IS  SUBPLOT  #  1$',  100,  3*.  D 


CALL  ENOGR  (0) 

CALL  OREL  (10.,  0.) 
CALL  AREA2D  (8.,  6.) 


{  Level  2  to  Level  1  } 


{  Level  1  to  Level  2  } 


CALL  FRAME 


CALL  HEADIN  ('THIS  IS  SUBPLOT  #  2$',  100,  3-.  D 


CALL  ENDGR  (0) 
CALL  DONEPL 


{  Level  2  to  Level  1  } 
{  Level  1  to  Level  0  } 


NOTE :  Programs  are  for  VAXcluster  only! 


■80 
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Example  2 


y.»*:y,y 
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Two  Dimensional  Graph i nc 


Once  the  subplot  area  has  been  defined,  the  next  step  is  to 
define  the  axis  system.  Below  are  some  of  the  major  routines, 
refer  to  the  User's  Manual  for  a  more  complete  description. 


Subroutine 


Level 

at 

call 


Level 

after 

call 


Primary  graph  setup  routines: 

CALL  GRAF  (XORG.  XSTP,  XMX,  YORG,  YSTP,  YMX) 
CALL  XLOG  (XORIG,  XCYCLE,  YORIG,  YSTEP) 

CALL  YLOG  (XORIG.  XSTEP,  YORIG.  YSTEP) 


CALL 

LOGLOG  (XORIG, 

XCYCLE,  YORIG.  YCYCLE) 

2 

3 

CALL 

POLAR  (THE  FAC, 

RSTEP,  XO 1  ST,  YD  1  ST) 

2 

3 

CALL 

GRID  (IXGRIO. 

IYGRID) 

3  * 

Same 

position  messages  in  INCHES: 

CALL 

MESSAG  (LMESS , 

IMESS,  XPOS,  YPOS) 

2.3 

Same 

CALL 

REALNO  (ANUM, 

1  PLACE,  XPOS.  YPOS) 

2.3 

Same 

CALL 

INTNO  (INUM,  XPOS,  YPOS) 

2.3 

Same 

CALL 

VECTOR  (XFROM, 

YFROM,  XTQ ,  YTO ,  1 VEC) 

2.3 

Same 

position  message  in  DATA  units: 

CALL 

RLMESS  (LMESS. 

IMESS,  XVAL,  YVAL) 

3 

Same 

CALL 

RLREAL  (ANUM. 

1  PL  ACE ,  XVAL.  YVAL) 

3 

Same 

CALL 

RLINT  (INUM.  XVAL ,  YVAL) 

3 

Same 

CALL 

RLVEC  (XFROM. 

YFROM.  XTO,  YTO,  IVEC) 

3 

Same 
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Now  let  us  reconsider  example  #2  as  it  actually  appeared  when 
ran.  Example  #3  is  the  same  program  with  the  sequence  of  commands 
used  to  draw  the  vectors  (with  arrowheads)  and  *he  integer 
distances  included.  It  is  actually  three  plots  (the  first  two  are 
thi  subplots,  the  third  is  the  entire  page)  so  that  the  vectors 
can  be  placed  around  the  subplot  areas  but  not  off  of  the  page. 
First  the  vector  is  drawn  (in  inches  from  the  origin),  then  INTNO 
places  the  value  currently  into  the  variable  INUM  (integer)  at 
the  coordinates  specified  (in  inches  from  the  origin).  Then 
MESSAG  is  called  to  place  the  inch  symbol  (notice  that  'ABUT1  is 
used  for  both  coordinates)  directly  after  the  previous  symbol  it 
just  wrote,  in  this  case  directly  after  what  INTNO  wrote. 


■x&m 
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PROGRAM  EXAMPLE 3 
CALL  COMPRS 
CALL  BGNPL  (0) 

CALL  PAGE  (24.,  11.) 

CALL  PHYSOR  (3- ,  2.) 

CALL  AREA2D  (8.,  6.) 

CALL  FRAME 

CALL  HEAD  I N ( ' TH I S  IS  SUBPLOT  #  1$‘ 
.100,  3-.  1.) 

CALL  ENDGR  (0) 

CALL  OREL  (10.,  0.) 

CALL  AREA2D  (8.,  6.) 

CALL  FRAME 

CALL  HEAOIN  ('THIS  IS  SUBPLOT  #  2$' 
,100,  3-.  1.) 

CALL  ENDGR  (0) 

CALL  PHYSOR  (0.,  0.) 

CALL  AREA2D  (20.,  11.) 

CALL  HEIGHT  (.5) 

CALL  SETCLR  ('RED') 

CALL  VECTOR  (3-.  0.,  3-.  2.,  3402) 

I NUM-2 

CALL  INTNO  (INUM,  3-5.  -75) 

CALL  MESSAG  (""  ,2, ' ABUT' , ' ABUT*) 
CALL  VECTOR  (0.,  4.,  3.,  4.,  3402) 

I NUM-3 

CALL  INTNO  (INUM,  1 .5,  4.5) 

CALL  MESSAG  (’ " ' , 2 , ' ABUT ' , 1  ABUT' ) 
CALL  VECTOR  (3-,  8.,  3-.  11., 3402) 
INUM-3 

CALL  INTNO  (INUM,  1 .5,  9.5) 

CALL  MESSAG  ( ' " 1 , 2 , ' ABUT 1 , 1  ABUT ' ) 
CALL  VECTOR  (3.,  3-,  11.,  3-, 3402) 

I NUM-8 

CALL  INTNO  (INUM,  6.5,  3-5) 

CALL  MESSAG  (""  ,  2,  'ABUT1,  'ABUT 


CALL  VECTOR  (4. , 2 . ,4. ,8. , 3402) 
INUM-6 

CALL  INTN0ONUM.4.5.5.5) 

CALL  MESSAG ( ""  , 2 , ' ABUT ' , ' ABUT ' ) 
CALL  VECTOR (11., 4., 13., 4., 3402) 

I NUM-2 

CALL  INTN0(INUM,11.75,4.5) 

CALL  MESSAG (""  ,2, 'ABUT' , 'ABUT') 

CALL  VECTOR (14. ,2. ,14. ,8. ,3402) 
INUM-6 

CALL  INTNO  (INUM, 14.5,5.5) 

CALL  MESSAG  (""  ,2, 'ABUT' , 'ABUT') 
CALL  VECTOR(13.,3.,21.,3.,3402) 

I NUM-8 

CALL  INTNO(INUM, 14.5,5.5) 

CALL  MESSAG  (""  ,2, 'ABUT' , 'ABUT') 
CALL  VECTOR  (21., 0., 21., 2., 3402) 

I NUM-2 

CALL  INTN0(INUM,21.5,.75) 

CALL  MESSAG (""  ,2, 'ABUT* , 'ABUT') 
CALL  VECTOR  (21  ., 8., 21  ., 11., 3402) 
INUM-3 

CALL  INTN0(INUM,21 .5,9.5) 

CALL  MESSAG (""  ,2, 'ABUT' , 'ABUT') 
CALL  VECTOR (21 . , 4 . , 24 . , 4 . , 3402) 
INUM-3 

CALL  INTN0(INUM,22.5,4.5) 

CALL  MESSAG (""  ,2, 'ABUT' , 'ABUT') 
CALL  ENDPL(O) 

CALL  DONEPL 
END 


) 
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TH8  18  SUBPLOT  #  1 


TH8 18  8UBPL0T  #  2 
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Once  OlSSPLA  is  in  level  3,  plotting  of  curves  may  begin. 
The  call  to  CURVE  may  be  made  as  many  times  as  needed. 


Subrout i ne 


Level  Level 
at  after 
cal  I  call 


To  draw  a  curve: 

CALL  CURVE  (XARAY,  YARAY,  NPTS.  I  MARK)  3  Same 

XARAY  -  array  of  X  values 
YARAY  -  array  of  Y  values 
NPTS  -  number  of  points  to  plot 

1  MARK  -  <  0  symbols  every  1th  point,  not  connected 
*  0  points  connected,  no  symbols 
>  0  symbols  every  I th  mark,  connected 

To  blank  an  area  around  the  curve: 


CALL  BLCURV  (OFWDTH,  OFLNTH) 

To  thicken  curves: 

CALL  THKCRV  (THKNSS) 

To  change  the  grace  margin  of  subplot: 

CALL  GRACE  (GRACEM) 

To  select  a  marker: 

CALL  MARKER  (ISYM) 

To  blank  symbol  markers: 

CALL  BLSYM 

To  change  size  of  markers: 

CALL  SCLPIC  (FACTOR) 


1,2,3  p/s 

1.2.3  P/s 

1.2.3  P/s 

1.2.3  P/s 

1.2.3  P/s 

1,2,3  p/s 
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This  example  reads  data  points  into  two  arrays,  then  plots 
the  curve.  This  object  must  be  drawn  to  scale,  so  PAGE  is  set  to 
23"  X  11".  AREA20  then  sets  the  plot  si2e  to  20"  X  8"  (from 
within  TITLE).  SETDEV  directs  error  messages  and  summary 
messages  to  unit  #  3  (F0R0O3.DAT  on  the  VAXcluster;  TAPE#  on  the 
CYBERs) ,  so  that  they  will  not  appear  on  the  screen  and  a 
hardcopy  is  kept.  INTAXS  will  print  only  the  integer  values  along 
the  axes,  and  YAXANG  will  write  the  values  along  the  Y  axis 
horizontally.  GRAF  sets  up  the  scale  of  the  axes,  the  X  axis  will 
run  from  0  to  20  by  steps  of  1,  the  Y  axis  will  run  from  -A  to  +4 
by  steps  of  1.  CURVE  is  then  called  twice,  once  to  plot  the  upper 
portion,  and  the  second  time  to  plot  the  lower  portion. 
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PROGRAM  EXAMPLES 

DIMENSION  X  (10000) ,Y  (10000) 

OPEN  (1 ,F ILE-' F0R001 .DAT' , STATUS-' OLD ' ) 
OPEN  (2 , F I LE* ' F0R002.DAT' , STATUS- ' OLD ' ) 


{  Level  0  to  Level  1  } 


CALL  COMPRS 
CALL  BGNPL  (0) 

CALL  SETDEV  (3,  3) 

CALL  SETCLR  ('BLUE') 

CALL  PAGE  (23-0.  11.0) 

CALL  HEIGHT  (.3) 

CALL  SWISSB 

CALL  SHDCHR  (90 . ,  1,  .002,  1) 

CALL  INTAXS 
CALL  YAXANG  (0.) 

CALL  TITLE  ( 1 MLTA  FOIL  PROFILES',  100,  'X  (  INCHES  )$', 
100,  'Z  (  INCHES  )$',  100,  20.,  8.)  {  Level  1  to  Level  2  } 

CALL  GRAF  (0.,  1.,  20.,  -U. ,  1.,  L.) 

{  Level  2  to  Level  3  ) 

CALL  SETCLR  ('RED') 

DO  20  1-1 , 10000 

READ(1,  '  (E  15*5. 3X.E  15*5)  '  .  END-21)  X  (I)  ,Y(I) 

CONTINUE 


CALL  CURVE  (X,  Y,  95OO,  0) 

DO  30  1-1 , 10000 

READ  (2,'  (E 15.5, 3X,F 15.5)  ' . END-32) X  ( I ) , Y ( I ) 


CONT 

INUE 

CALL 

CURVE 

(X,  Y,  9500,  0) 

CALL 

ENDPL 

(0) 

{  Level 

3  to  Level 

1  } 

CALL 

DONEPL 

{  Level 

1  to  Level 

0  } 

END 
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Legends  and  Stor i es 

.  To  draw  a  Legend  for  your  graph,  or  write  out  a  block  of 
text  (Story),  use  one  of  the  following  routines. 


Subrout i ne 

Level 

Level 

at 

after 

cal  1 

cal  1 

To 

store  1 i nes  of  text: 

■ 

CALL  LINES  (LSTRING,  IPKARAY,  1 L 1 NE) 

1.2,3 

Same 

To 

draw  a  legend: 

CALL  LEGEND  (IPKARAY,  NUNES,  XPOS, 

YPOS) 

1,2,3 

Same 

To 

initialize  IPKARAY  and  set  J  to  max 

1  i  nes 

possible 

i n  Legend 

J  -  LINEST  (IPKRAY,  NPKWRD,  IMAX) 

1.2,3 

Same 

To 

write  a  block  of  centered  text: 

CALL  STORY  (IPKARAY,  NLINES,  XPOS,  YPOS) 

1.2,3 

Same 

To 

left  justify  a  block  of  text: 

CALL  LSTORY  (IPKARAY,  NLINES,  XPOS, 

YPOS) 

1,2,3 

Same 

To 

right  justify  a  block  of  text: 

CALL  RSTORY  (IPKARAY,  NLINES,  XPOS, 

YPOS) 

1.2,3 

Same 

To 

get  the  size  of  LEGEND: 

XLEN  -  XLEGEND  (IPK,  NLINES) 

1.2,3 

Same 

YLEN  -  YLEGEND  (IPK,  NLINES) 

1.2,3 

Same 

To 

get  the  size  of  STORY: 

XLEN  =  XSTORY  (IPK,  NLINES) 

1,2,3 

Same 

YLEN  =  YSTORY  (IPK,  NLINES) 

1.2,3 

Same 
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To  get  the  length  of  messages: 

XLEN  -  XMESS  (LMESS,  I  MESS)  1,2,3 

To  get  the  length  of  numbers: 

XLEN  -  XREAL  (ANUM,  IPLACE)  {  REALS  }  1,2,3 

XLEN  -  XI NT  (I NT)  {INTEGER}  1,2,3 

To  get  X,Y  position  in  inches  of  a  data  point: 

A  -  XPOSN  (XVAL,  YVAL)  1,2,3 

B  -  YPOSN  (XVAL,  YVAL)  1,2,3 

To  get  X, Y  coordinate  values: 

A  *  XINVRS  (XINCH,  YINCH)  1,2,3 

B  -  YINVRS  (XINCH,  YINCH)  1,2,3 


Same 


Same 

Same 


Same 

Same 


Same 

Same 


The  next  example  uses  a  legend  containing  9  lines  which  are 
packed  by  LINES  and  stored.  MYLEGN  changes  the  title  on  the 
legend  to  any  string.  BLREC  uses  the  values  of  XLEGEND  and 
YLEGEND  to  blank  out  the  area  where  the  legend  will  be  written, 
to  prevent  curves  from  over  writing  it.  Notice  that  LEGLIN  is 
called  directly  before  CURVE,  the  line  type  is  stored  with  the 
legend  titles.  NOTE :  the  curves  must  be  called  in  the  same  order 
as  they  appear  in  the  legend.  LEGEND  then  prints  out  the  legend 
at  the  specified  coordinates. 

The  second  example  is  basically  the  same  except  SPCMOD  is 
called.  This  subroutine  calls  MYSPEC  and  also  stores  the  color 
of  the  line  that  was  used. 
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PROGRAM  EXAMPLES 

DIMENSION  IPAK (10000) ,Y (10) ,2(10) 
OPEN  (1 , F I LE- ' JETDATA ' , STATUS- ' OLD ' ) 


CALL 

COMPRS 

CALL 

BGNPL 

(0) 

CALL 

SETDEV 

(2.  2) 

CALL 

SETCLR 

('BLUE') 

CALL 

PAGE 

(12.0,  11.0) 

CALL 

SWISSB 

CALL 

INTAXS 

CALL 

YAXANG 

(0.) 

CALL 

SHDCHR 

(90..  1,  .002,  1) 

CALL 

TITLE 

(' IMPACT  WINDOWS' , 

100,  'Z  - 

.  *  Y  -  1 

IMPACTS' 

,  100,  9..  7.) 

CALL 

GRAF 

(-200.,  20.,  200., 

-200.,  20 

CALL 

SETCLR 

('RED') 

CALL 

LINES 

(' .0000  -  .0625$'  . 

IPAK,  1) 

CALL 

LINES 

(' .0625  "  .1220$' , 

IPAK,  2) 

CALL 

LINES 

('.1250  -  .2500$', 

IPAK,  3) 

CALL 

LINES 

(' .2500  -  .5000$' , 

IPAK,  4) 

CALL 

LINES 

('.5000  -  .7500$', 

IPAK,  5) 

CALL 

LINES 

C .7500  -  1.000$' , 

IPAK,  6) 

CALL 

LINES 

(' 1 .000  -  2.000$' , 

IPAK,  7) 

CALL 

LINES 

(' 2.000  -  3.000$', 

IPAK,  8) 

CALL 

LINES 

(' 3.000  -  A. 000$ ' , 

IPAK,  9) 

CALL 

LINES 

('  >  4.000$',  IPAK, 

,  10) 

CALL 

MYLEGN 

('DEPTH  IN  INCHESS 

',  100) 

XR  - 

XLEGND 

(IPAK,  10)  +  .3 

YR  - 

YLEGND 

(IPAK,  10)  +  .3 

CALL 

BLREC 

(7. -.3,  6. -.3,  XR+ 

.3,  YR+. 3, 

DO  20  1-1,500 

MPACT$ 

200.) 


100, 


.02) 


READ  (1 . '  (70X.F5.2, 
IF  ((DEP.GE.0.0  ) 

IF  ((DEP.GE. 0.0625) 
IF  ((OEP.GE. 0.1250) 
IF  ((DEP.GE. 0.25  ) 
IF  ((OEP.GE. 0.5 
IF  ((DEP.GE. 0.75 
IF  ((DEP.GE. 1 .0 
IF  ((DEP.GE. 2.0 
IF  ((DEP.GE. 3-0 
IF  ((DEP.GE. L.O) 
CALL  LEGLIN 
CALL  CURVE  (Y, 
CONTINUE 
CALL  RESET 
LEGEND 
ENDPL 
DONEPL 


.AND. 

.AND. 

.AND. 

.AND. 

.AND. 

.AND. 

.AND. 

.AND. 

.AND. 

CALL 


(DEP.LT. 0.0625) 
(DEP.LT. 0.1250) 
(DEP.LT. 0.2500) 
(DEP.LT. 0.5000) 
(DEP.LT. 0.750  ) 
(DEP.LT. 1 .0 
(DEP.LT. 2.0 
(DEP.LT. 3.0 
(DEP.LT. L.O 
MARKER  (10) 


,  END-32) DEP 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 


) 

) 

) 

) 

) 

)) 

)) 

)) 

)) 


,  Y  (1)  ,Z  (1) 
MARKER  (1) 
MARKER (2) 
MARKER  (3) 
MARKER  (4) 
MARKER (5) 
MARKER  (6) 
MARKER  (7) 
MARKER  (8) 
MARKER  (5) 


( 1 BLNKS ' ) 
(IPAK,  10, 
(0) 


tit 
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PROGRAM  EXAMPLE6 

DIMENSION  X  (50)  ,Y  (50)  .IPAK  (1000) 

OPEN  ( 1 , F I LE- '  DAT 1 , ST ATUS- 1  OLD 1 ) 

CALL  COMPRS 
CALL  BGNPL  (0) 

CALL  SETCLR  ('BLUE') 

CALL  HEIGHT  (.2) 

CALL  PAGE  (12.0,  11.0) 

CALL  INTAXS 
CALL  YAXANG  (0.) 

CALL  TITLE  ('FUNCTIONS',  100,  'X  -  AXISS',  100, 
.'Y  -  AXISS',  100.  9.,  7-) 

CALL  SETCLR  ('GREEN') 

CALL  GRAF  (-3.,  1.,  3*.  -*»•.  1  •  •  10.) 

CALL  SETCLR  ('GREEN') 

CALL  LINES  ('INPUT  POINTSS',  IPAK,  1) 

CALL  LINES  ('  OEGREE  1  $'.  IPAK,  2) 

CALL  LINES  ('  DEGREE  2  $',  IPAK,  3) 

CALL  LINES  ('  DEGREE  3  S',  IPAK.  A) 

CALL  LINES  ('  DEGREE  A  $'.  IPAK,  5) 

CALL  LINES  ('  DEGREE  5  S',  IPAK.  6) 

CALL  LINES  ('  DEGREE  6  S'.  IPAK.  7) 

CALL  MYLEGN  ('POLY  DEGREES'.  100) 

XR  «  XLEGND  (IPAK,  7)  +  -3 
YR  -  YLEGND  (IPAK,  7)  +  .3 
CALL  BLREC  (7.-. 3,  6.-. 3,  XR+.3,  YR+.3,  .02) 

DO  10  1-1,13 
READ  (1  ,*)  X  (I)  ,  Y  (I) 

10  CONTINUE 

CALL  SPCMOD 
CALL  LEGLIN 

CALL  CURVE  (X,  Y,  13,  1) 

DO  30  J-1,6 
DO  20  1-1, A9 

READ  (1  ,*)  X  (I)  ,  Y  (I ) 

20  CONTINUE 

CALL  SPCMOD 
CALL  LEGLIN 

CALL  CURVE  (X,  Y,  A9,  1) 

30  CONTINUE 

32  CALL  BLMOVE  (XR+.9,  0.) 

CALL  LEGEND  (IPAK,  7,  7-,  6.) 

CALL  ENDPL  (0) 

CALL  DONEPL 
END 

SUBROUTINE  MYSPEC  (J) 

IF  (J  .EQ.  1  )  CALL  SETCLR  ('RED') 

IF  (J  .EQ.  2  )  CALL  SETCLR  ('MAGENTA') 

IF  (J  .EQ.  3  )  CALL  SETCLR  ('BLUE') 

IF  (J  .EQ.  I*  )  CALL  SETCLR  ('CYAN') 

IF  (J  .EQ.  5  )  CALL  SETCLR  ('YELLOW') 

RETURN 

END 


i 

-> 

"1 
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InterDotat ion 

To  use  interpolation  rout 

ines  provided  by 

DISSPLA: 

Subrout i ne 

Level 

Level 

CALL  SPLINE 

at 

cal  1 

1.2,3 

after 
cal  1 

P/s 

CALL  PSPLIN 

1.2.3 

p/s 

CALL  P0LY3 

1.2.3 

P/s 

CALL  P0LY5 

1.2.3 

P/s 

CALL  LINEAR 

1.2.3 

P/s 

CALL  PARA3 

1.2,3 

P/s 

CALL  PARA5 

1.2.3 

P/s 

CALL  STEP 

1.2,3 

p/s 

CALL  BARS  (BARWTH) 

1.2,3 

P/s 

For  smoothing:  (Supply  error  weighting  factors 

in  blank 

common  YDLARAY) 

CALL  SMOOTH 

1.2,3 

P/s 

CALL  PSMTH 

1,2,3 

p/s 

To  change  line  texture: 

CALL  DOT 

1.2,3 

P/s 

CALL  DASH 

1,2,3 

P/s 

CALL  CHNDOT 

1,2,3 

P/s 

CALL  CHNDSH 

1,2,3 

P/s 

CALL  RESET  ('DOT')  {SOLID} 

1,2,3 

P/s 

To  construct  your  own  texture: 

CALL  MRSCOD  (TLENG,  NMRKSP, 

\ 

RATRAY) 

1.2,3 

P/s 
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PROGRAM  EXAMPLE 7 

DIMENSION  X  ( 1 3)  , Y  ( 1 3)  , RATRAY  (6) 

DATA  X/-3.. -2. 5. -2.. -1. 5.-1.. -.5.0... 5.1.. 1.5. 2., 
Y/" 1-»i**5.-3..6..-2.,3.,-3.,8.,-2.,5..-i.,5. 


DATA  RATRAY  /5.,  2. 
CALL  COMPRS 
CALL  BGNPL  (0) 
CALL  HEIGHT  (1.) 
CALL  SCMPLX 
CALL  PAGE 
DO  30  J-1,6 
CALL  PHYSOR 
CALL  ARE A2D 
CALL  SETCLR 
CALL  GRAF 
CALL  FRAME 
CALL  MESSAG  (' 


2..  2..  7*.  2./ 


(24.0,  18.0) 

(1..  FLOAT  (J-l)*3.) 

(15*.  3-) 

(‘GREEN1) 

(-3..  1..  3*.  -4..  ; 


10.) 


IF  (  J.EQ.l  )  THEN 
CALL  LINEAR 
CALL  SETCLR  ('RED') 
CALL  MESSAG  ('Ll NEAR$ 1 

ELSE  IF  (  J.EQ.2  )  THEN 
CALL  P0LY3 
CALL  SETCLR  ('BLUE') 
CALL  MESSAG  ( ' P0LY3$ ‘ . 
CALL  OOT 

ELSE  IF  (  J.EQ.3  )  THEN 
CALL  STEP 
CALL  SETCLR 
CALL  MESSAG 
CALL  DASH 

ELSE  IF  (  J.EQ.4  )  THEN 
CALL  SPLINE 


100.  15-.  0.19) 


100,  16.,  1.) 


100,  16.,  1.) 


('YELLOW') 

('STEPS’,  100,  16.,  1.) 


2.5,3./ 

.-3./ 


CALL  SETCLR  ('MAGENTA') 

CALL  MESSAG  ('SPLINES',  100,  16.,  1.) 
CALL  CHNOOT 

ELSE  IF  (  J.EQ.5  )  THEN 
CALL  BARS  (1.) 

CALL  SETCLR  ('CYAN') 

CALL  MESSAG  ( ' BARS$ ' ,  100,  16.,  1.) 
CALL  CHNDSH 

ELSE  IF  (  J.EQ.6  )  THEN 
CALL  P0LY5 

CALL  SETCLR  ('YELLOW') 

CALL  MESSAG  ('P0LY5S ' ,  100,  16.,  1.) 
CALL  MRSCOD  (2.,  6,  RATRAY) 

END  IF 

CALL  CURVE  (X,  Y,  13,  1) 

CALL  ENOGR  (0) 

CONTINUE 

CALL  ENOPL  (0) 

CALL  DONEPL 
END 


0I5SPLA 


P3- 


P0LY5 

BARS 


SPLINE 


P0LY3 


LINEAR 
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Example  #8  shows  a  polar  plot  of  flows  around  a  prop.  The 
data  is  read  in  from  two  data  files  containing  angles  and 
magnitudes  of  the  flows.  A  vector  will  be  drawn  for  each  flow, 
will  the  length  being  the  magnitude,  the  direction  being  the 
angle.  The  data  did  have  to  be  adjusted  to  convert  to  the  polar 
coordinate  system.  Notice  that  subroutine  TITLE  is  used  even 
though  it  is  is  an  obsolete  routine  it  is  still  supported  by 
DISSPLA.  In  ail  calls  to  MESSAG,  the  position  (XVAL,  YVAL)  is 
given  in  polar  coordinates  (i.e.,  ANGLE,  MAGNITUDE).  In  the  call 
to  POLAR  you  determine  whether  angles  are  in  radians  or  degrees. 


m 
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PROGRAM  EXAMPLE8 

OPEN  (1 ,F ILE-'TUBEI .DAT’ .STATUS-' OLD') 

OPEN  (2 , F I LE- ' TUBE2 . DAT  * , STATUS- ' OLD ' ) 

CALL  COMPRS 
CALL  SETDEV  (9,  9) 

CALL  BGNPL  (1) 

CALL  PAGE  (U.,  II.) 

CALL  SETCLR  ('BLUE') 

CALL  TITLE  ('  ', 1 , 'RADIUS' ,6, 1H  ,1,10.5,10.5) 
CALL  POLAR  (3- 14159/180. , .25,5.25,5.25) 

CALL  HEIGHT  (.2) 

CALL  SETCLR  ('GREEN') 

CALL  GRID  (1,1) 

DO  105.K-1.2 
IF(K.EQ.I)  THEN 
CALL  SETCLR  ('RED') 

READ  (1 ,  1  (5X,  1 8X .  13)  ')  INUM 

CALL  RLMESS  ('  TUBE  -  ',8,60., 1.3) 

CALL  INTNO  ( I NUM, ' ABUT 1 , ' ABUT ' ) 

CALL  RLMESS  ('  TUBE  -  ' , 8, 1 35 . . 1 . 6) 

CALL  I NTNO  ( I NUM, ' ABUT ' , ' ABUT ' ) 

CALL  RLMESS  ('  TUBE  -  ' ,8,220. , 1 .6) 

CALL  INTNO  ( I NUM, ' ABUT ' , ' ABUT 1 ) 

CALL  RLMESS  ('  TUBE  -  ' . 8. 310 . , 1 . 35) 

CALL  I NTNO  ( I NUM, 1  ABUT ' , ' ABUT ' )  50 

endif  105 

IF  (K . EQ . 2)  THEN  999 

CALL  SETCLR  ('BLUE ') 

READ  (2,'  (5X.18X.I3)  ')  INUM 

CALL  RLMESS  ('  TUBE  -  ‘,8,50. ,1.3) 

CALL  INTNO  ( I NUM, ' ABUT ' , ' ABUT ' ) 

CALL  RLMESS  ('  TUBE  -  1 ,8. 140. , 1 .5) 

CALL  INTNO  ( I NUM, ' ABUT ' , 1  ABUT  1 ) 

CALL  RLMESS  ('  TUBE  -  ' , 8, 220 . , 1 . 8) 

CALL  INTNO  ( I NUM, 1  ABUT ' , ' ABUT ' ) 

CALL  RLMESS  ('  TUBE  -  ' , 8 , 290 . , 1 . 3) 

CALL  INTNO  ( I NUM, ' ABUT ' , ' ABUT ' ) 

ENDIF 


I F  ( (TH . GT . 0)  .AND.  (TH.LE.90))  THEN 
IF  (THNEW.LT. 0)  THNEW— THNEW 
ENDIF 

IF  ((TH.GT.180) .AND.  (TH .LE . 270) ) THEN 
ANG  -  2. 

XI  -  R1*C0SD  (TH-TH1-ANG* (90.) ) 

Y1  -  R1*SIND  (TH-TH1-ANG* (90.)) 

IF  ( (TH.GT . 180)  .AND.  (TH . LE . 270) ) THEN 
IF  ((Xl.LT.O.)  .AND.  (Yl.LT.O.))  THEN 
Y 1  -  -Y  1 

XI  -  -XI 

ENDIF 

XNEW  -  X  -  XI 
YNEW  -  Y  -  Y1 

RNEW  -  SQRT((XNEW**2)  +  (YNEW**2)) 
THNEW  -  ATAND  (YNEW/XNEW) 

IF  (TH.GT. 180)  .AND.  (TH.LE270) )  THEN 
IF  (THNEW. GT.O)  THNEW  -  THNEW+180. 
ENDIF 
ENDIF 

CALL  RLVEC (TH,R, THNEW, RNEW, 1201) 

CONTINUE 

CONTINUE 

CONTINUE 

CALL  ENDPL(O) 

CLOSE  (9. STATUS- 'DELETE ' ) 

CALL  DONEPL 
END 


00  50  1-1,2500 

IF  (K.EQ. 1)  READ  (1 , 1 1 . END-105) TH , R 1 ,TH 1 
IF  (K.EQ.2)  READ  (2,1 1 , END-999) TH , R 1 ,TH 1 
IF  (K.EQ. 1)  R  -  1.25 
IF  (K.EQ.2)  R  -  .75 

X  =  R  *  COSD(TH) 

Y  =  R  *  S I  NO  (TH) 

XI  *  R1  *  COSO  (180.  +TH  -  TH 1 ) 

Y 1  =  R1  *  SIND  (180.  +  TH  -  TH1) 

XNEW  -  X  +  XI 

YNEW  -  Y  +  Y 1 

RNEW  -  SQRT((XNEW**2)  +  (YNEW**2)) 

THNEW  -  ATAND (YNEW/XNEW) 
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Post-Processor 

When  COMPRS  is  used  as  the  device,  DISSPLA  creates  a  device 
independent  META  file  (called  DISPLOT.DAT  on  VAX  and  PLFILE  on 
the  CYBERS) .  When  you  run  the  post  processor,  you  can  choose  any 
device  and  make  certain  modifications.  A  plot  can  be  scaled  up  or 
down,  certain  ones  drawn,  or  a  small  portion  of  the  plot  "zoomed" 
in  on. 

The  following  command  used  on  the  preceding  example  would  produce 
the  following  graph: 


MODI  -  1  (  WINDOW  -  UPPER  (  5-75). LOWER  (  0.75), 

RIGHT(  5. 75). LEFT  (  0.75)  *  SIZE  -  11,8) 


NOTE:  The  graph  cannot  be  changed  (i.e.,  viewpoint),  you  can  only 
work  with  the  plot  as  it  exists! 
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3-Diwensional  Plotting 


Three-dimensional  plotting  in  DISSPLA  is  similar  to  two- 
dimensional  plotting  except  that  some  subroutines  are  changed  to 
accommodate  the  Z  coordinate  arguments.  There  are  alto  some 
additional  routines  to  cope  with  3~D  plotting  properties.  The 
user  wishing  to  use  3“D  is  advised  to  read  the  User's  Manual. 


Subrout i ne  Level  Level 

at  after 

cal  1  cal  1 

To  select  relative  axis  lengths: 

CALL  V0LM3D  (X  3 AX  I S ,  Y3AXIS.  23AXIS)  2  p/s 

To  select  a  viewpoint:  (Absolute  units) 

CALL  VUABS  (XVU,  YVU,  ZVU)  2  p/s 

CALL  VUANG  (PHI.  THETA,  RADIUS)  2  p/s 

To  set  3~D  axis: 

CALL  GRAF  3D  (X30RIG,  X3STEP,  X3MAX, 

Y30RIG,  Y3STEP,  Y 3MAX ,  2  3 

Z30RIG,  Z3STEP,  Z3MAX) 

To  draw  a  curve: 

CALL  CURV3D  (XARAY,  YARAY,  ZARAY,  N.  I)  3  Same 

N  -  Number  of  points  to  be  plotted 

I  -  <  0  Symbol  every  i-th  mark  not  connected 
■0  No  symbols  just  a  line 
>  0  Symbol  every  i-th  msrk  and  connected 


To  draw  3~D  Vector: 

CALL  VECTR3  (XFROM,  YFROM,  ZFROM, 

XTO,  YTO,  ZT0,  I VEC)  2,3  Same 


CALL  RLVEC3  (XROM,  YFROM,  ZFROM, 
XTO,  YTO,  ZTO,  I VEC) 


2.3 


Same 
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Example  #10  shows  a  simple  3-D  plot  of  a  conic  section. 
Notice  that  along  with  AREA2D,  AXES3D  is  called  to  define  the  3“& 
work  area.  It  is  defined  in  "absolute  units",  so  that  VUABS  can 
define  the  view  angle  relative  to  the  work  area.  Four  subplots 
will  be  put  onto  the  same  page,  each  one  using  a  different 
viewpoint.  The  X,  Y,  and  Z  arrays  are  filled  with  the  coordinates 
to  be  plotted  using  the  equation  of  the  conic. 
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PROGRAM  EXAMPLE  10 

DIMENSION  X  (500) ,  Y (500)  ,  Z  (500) ,  ZX 
CALL  COMPRS 
CALL  BGNPL  (0) 

CALL  PAGE  (14.,  11.) 

CALL  HEIGHT  (.3) 

CALL  INTAXS 
CALL  ZAXANG  (90.) 

DO  100  I  RUNS* 1 ,4 
CALL  SETCLR  ('BLUE') 

IF  (IRONS. EQ.l)  THEN 
CALL  PHYSOR  (.5,  .5) 

CALL  AREA2D  (6.,  5.) 

CALL  VUABS  (3..  -30.,  3-) 

ELSE  IF  (I RUNS. EQ. 2)  THEN 
CALL  PHYSOR  (.5,  5-5) 

CALL  AREA2D  (6.,  5.) 

CALL  VUABS  (1.5,  1 .5,  30.) 

ELSE  IF  (IRUNS.EQ.3)  THEN 
CALL  PHYSOR  (7-5,  -5) 

CALL  AREA2D  (6.,  5.) 

CALL  VUABS  (30.,  30.,  30.) 

ELSE  IF  ( I  RUNS . EQ.4)  THEN 
CALL  PHYSOR  (7-5.  5*5) 

CALL  AREA2D  (6.,  5.) 

CALL  VUABS  (-30.,  -30.,  *30.) 

END  I  F 

CALL  AXES3D  {'X1,  1,  ' Y«  ,  1,  'Z\  1,  3-.  3-.  3.) 

CALL  GRAF3D  (-30.,  15.,  30.,  -30.,  15- .  30.,  -30.,  15.,  30.) 
DO  10  1-1,60 
I  COUNT  -1 
DO  20  J-1,60 
ZZ  -  FLOAT  (1-30) 

YY  -  FLOAT  (J-30) 

XX  -  (((ZZ**2/l6)-(YY**2/4))  +  1  )  *4 
IF(XX.GE.O)  THEN 
X  (I  COUNT)  -  SQRT(XX) 

Y  (I  COUNT)  -  YY 
Z  (I  COUNT)  -  ZZ 
ZX  ( I  COUNT)  -  -SQRT(XX) 

I  COUNT  -  I COUNT+1 
END  IF 

20  CONTINUE 

I  COUNT  -  I  COUNT- 1 
CALL  SETCLR  ('GREEN') 

DO  500  IX- 1,1  COUNT 
Z  ( I COUNT+ 1  X)  -  Z ( I COUNT+1 -  I  X) 

X  ( I COUNT+ 1  X)  -  ZX (I COUNT+l- IX) 

Y  (I  COUNT+1  X)  -  Y  (I  COUNT+1 -IX) 

CONTINUE 


Z (2* I COUNT+ 1 ) -Z  ( 1 ) 

X (2* I COUNT+ 1 ) -X  ( 1 ) 

Y  (2*1  COUNT+1) -Y(l) 
CALL  CURV3D  (X,  Y,  Z. 
.  2  *  I COUNT+1, 0) 

10  CONTINUE 
CALL  ENDGR  (0) 

100  CONTINUE 

CALL  ENDPL(O) 

CALL  DONEPL 
END 


500 
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Example  #11  shows  how  to  plot  a  surface  defined  by  a 
function  of  two  variables.  Once  the  work  box  is  set  up  VUANGL  is 
called  to  define  the  viewpoint.  This  has  the  same  effect  as 
VUABS.  only  it  uses  spherical  coordinates.  The  equation  to  be 
plotted  is  put  in  an  EXTERNAL  function  with  two  arguments.  The 
call  to  SURFUN  contains  the  name  of  the  function  and  will 
automatically  compute  the  points  to  be  plotted.  The  figure  is 
plotted  15  times,  each  from  a  different  viewpoint  to  simulate 
rotation. 


Example  #12  shows  a  single  surface  plot,  the  call  to  VUABS 
will  draw  the  surface  as  if  you  were  looking  at  it  from  (*• 
10. , U. ,20.)  .  ZAXANG  (0 .)  will  label  the  Z-axis  hori2ontally  the 
same  way  YAXANG  works  in  2-D. 


NOTE:  When  defining  the  viewpoint,  you  must  draw  the  object  from 
the  outside  looking  in.  DISSPLA  will  print  an  error  message  if 
you  are  within  the  work  box  (i.e.,  looking  from  the  inside  out). 
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This  page 


intentionally  left  blank. 
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PROGRAM  EXAMPLE  11 
EXTERNAL  PYRR,  PYRR1 
CALL  COMPRS 
CALL  BGNPL  (0) 

THETA  -  0 

PHI  -  -180. 

CALL  PAGE  (35-5.  30.) 

DO  20  1-1,5 
DO  30  J* 1,3 
CALL  SETCLR  ('BLUE') 

CALL  PHYSOR  (FLOAT  ( I - 1) *  (35 -5/5 •) .  FLOAT  (J) *10-10) 

CALL  AREA2D  (35-5/5 • .  30. /3.) 

CALL  FRAME 
CALL  BLSUR 

CALL  AXES3D  (0,  0,  0,  0,  0,  0,  5-.  5-.  5») 

CALL  VUANGL  (PHI,  THETA,  20.) 

CALL  GRAF  3D  (-1000.,  500.,  1000.,  -1000.,  500.,  1000.,  -1000., 

500., 1000.) 

CALL  SETCLR  ('CYAN') 

CALL  SURFUN  (PYRR,  1,  100.,  1,  100.,  WORK) 

CALL  SETCLR  (‘RED’) 

CALL  SURFUN  (PYRR1 ,  1,  100.,  1,  100.,  WORK) 

CALL  ENDGR  (0) 

IF  (  J.EQ.2)  THEN 
THETA  -  90.-  1*10. 

ELSE 

PHI  -  PHI  +  90. 

END  IF 
CONTINUE 
PHI  -  -l80. 

THETA  -  I  *  10. 

CONTINUE 
CALL  ENDPL  (0) 

CALL  OONEPL 
END 


FUNCTION  PYRR(X.Y) 

PYRR  -  SQRT(ABS  ( (3*Y**2-l6*X) /12) ) 

RETURN 

END 


FUNCTION  PYRR 1  (X.Y) 

PYRR 1  -  -SQRT(ABS((3*Y**2-16*X)/12)) 

RETURN 

END 


mm WPS? 
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PROGRAM  EXAMPLE  12 
EXTERNAL  PYRR 
CALL  COMPRS 
CALL  BGNPL  (0) 

CALL  PAGE  (11*.,  11.) 

CALL  HEIGHT  (.3) 

CALL  INTAXS 
CALL  ZAXANG  (90.) 

CALL  SETCLR  ('BLUE') 

CALL  AREA2D  (10..  8.) 

CALL  HEADIN  ('  3-D  SURFACE  PLOTS',  100,  1.2,  1) 
CALL  AXES3D  ('X'.  1,  ' Y ' ,  1.  'Z'.  1,  5*.  5-.  5-) 
CALL  VUABS  (-10..  1* . ,  20.) 

CALL  GRAF3D  (-2.,  1.,  2.,  -2.,  1.,  2.,  0.,  1.,  2.) 
CALL  SETCLR  ('RED') 

CALL  SURFUN  (PYRR,  2,  .1,  2,  .1,  WORK) 

CALL  ENDPL  (0) 

CALL  OONEPL 
END 

FUNCTION  PYRR(X.Y) 

PYRR  -  (X**2+2*Y**2) *EXP  ( 1 -X**2-Y**2) 

RETURN 

END 
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3~D  Surf  ace  Plot 
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Graff i t i  Plots 


Graffiti  plots  allow  the  user  to  project  a  2-D  plot  onto  a 
plane  in  3“0.  The  call  below  defines  the  plane  by  giving  three 
points  that  are  on  the  plane. 


CALL  GRFITI  (XLCORN,  YLCORN ,  ZLCORN, 
XBASEX,  YBASEX,  ZBASEX, 
XOTHER,  YOTHER.  ZOTHER) 


XLCORN,  YLCORN,  ZLCORN  -  coordinates  of  the  lower  left  corner  of 

the  plane  in  absolute  workbox  units. 

XBASEX,  YBASEX,  ZBASEX  -  coordinates  of  a  point  on  the  X-axis  in 

absolute  workbox  units. 

XOTHER,  YOTHER,  ZOTHER  -  coordinates  of  a  point  in  the  plane  above 

the  X-axis  in  absolute  workbox  units. 


Example  #13  shows  a  3~0  graffiti  plot.  It  is  important  to 
note  that  when  drawing  a  graffiti  plot,  the  2-D  plot  must  be 
finished  (  CALL  END3GR  )  before  3~0  plotting  or  another  2-D  plot 
can  be  started.  The  workbox  is  defined  to  be  5  X  5  X  5,  the  first 
graffiti  plot  defines  the  2-D  plane  (0,5,0)  (5,5,0)  (2,5,2)  which 
is  the  Z-X  plane  when  Y  is  equal  to  5*  A  grid  is  drawn  on  the 
plane  and  a  message  written  before  it  is  ended  by  END3GR.  Then 
the  second  plane  is  defined,  the  Y-Z  plane  when  X  is  5.  A  grid  is 
again  drawn  on  the  plane  and  the  message  written  before  it  is 
ended  by  END3GR.  Finaly  the  third  plane  is  defined,  the  X-Y  plane 
when  Z  is  zero,  and  the  same  events  occur.  After  the  final  call 
to  END3GR  the  3“D  plotting  can  begin. 


Example  #11*  shows  a  graffiti  plot  with  actual  graphs. 
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Th i s  page  intentionally  left  blank. 
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PROGRAM  EXAMPLE  13 


(0. 


CALL  TK41 
CALL  BGNPL 
CALL  PAGE 
CALL  SETCLR 
CALL  PHYSOR 
CALL  HEIGHT 
CALL  INTAXS 
CALL  AREA2D 
CALL  SETCLR 
CALL  FRAME 
CALL  AXES3D 
CALL  VUABS 
CALL  GRAF  3D 
CALL  GRFITI 
CALL  AREA20 
CALL  GRAF 
CALL  SETCLR 
CALL  MESSAG 
CALL  MESSAG 
CALL  BLREC 
CALL  GRID 
CALL  END3GR 
CALL  GRFITI  (5 
CALL  AREA2D 
CALL  GRAF 
CALL  SETCLR 
CALL  MESSAG 
CALL  MESSAG 
CALL  BLREC 
CALL  GRID 
CALL  END3GR 
CALL  GRFITI  (0 
CALL  AREA2D 
CALL  GRAF 
CALL  SETCLR 
CALL  MESSAG 
CALL  MESSAG 
CALL  BLREC 
CALL  GRID 
CALL  END3GR 
CALL  SETCLR 
CALL  BOX 30 
CALL  ENDPL 
CALL  DONEPL 
END 


(U105) 

(0) 

(1A.,  11.) 

('BLUE') 

(.5.  .5) 

(.25) 

(13..  10.) 

('RED') 

('X  -  AXIS' , 8 , * Y  -  AXIS' .8. *2  -  AXIS' ,8,5- .5- ,5-) 

(-10.,  -10.,  10.) 

(-150., 50., 150., -150., 50., 150.. -150., 50., 150.) 


150.,  -150., 

1., 


,  5- .  0. ,  ! 

(5.,  5.) 

(-150.,  10 

('GREEN') 

('THIS  IS  THE$ ' ,  100, 

('Z  -  X  PLANES',  100,  1. 
(.7,  1.7.  3..  1.5,  .08) 
(1,  1) 

(0) 

,  5*,  0.,  5*.  0.,  0.,  5. 
(5.,  5.) 

(-150., 10., 150. 

(‘YELLOW) 

('THIS  IS  THE$' ,  100, 

CY  -  Z  PLANES',  100, 


5 

10., 

2.5) 

2.) 


2.) 

150.) 


2.) 


150. ,10. ,150.) 


5) 

.) 


(.7.  1-7,  3.,  1.5,  .08) 

(1,  1) 

(0) 

,  0.,  0.,  5.,  0.,  0.,  0.,  5.,  0.) 

(5..  5.) 

(-150., 10., 150., -150., 10., 150.) 

('BLUE') 

('THIS  IS  THES ' ,  100,  1 .,  2.5) 

('X  -  Y  PLANES',  100,  1.,  2.) 

(.7.  1.7,  3..  1.5,  -08) 

(1,  1) 

(0) 

('CYAN') 

(0) 
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Character  Sets 

DISSPLA  contains  numerous 

alphabets 

and  character 

styles . 

Consult  the  User's  Manual  for 

examples 

of  each. 

Subrout i ne 

Level 

at 

cal  1 

Level 
after 
cal  1 

CALL  CARTOG 

1,2,3 

P/s 

CALL  SIMPLX 

1,2.3 

P/s 

CALL  SCMPLX 

1,2,3 

p/s 

CALL  COMPLX 

1,2,3 

p/s 

CALL  DUPLX 

1,2,3 

p/s 

CALL  TRIPLX 

1.2,3 

p/s 

CALL  GOTHIC 

1,2,3 

p/s 

Then  to  choose  an  alphabet: 

CALL  BASALF  (LALPHA) 

1,2,3 

p/s 

where  LALPHA  is: 

"STANDARD"  "L/CSTD"  "GREEK"  "L/CGREEK"  "RUSSIAN"  "L/CRUSSIAN"  "HEBREW" 
"ITALIC"  "L/CITALIC"  "SCRIPT"  "SPECIAL"  "MATHEMATIC"  "INSTRUCTION" 


Example  #15  shows  how  to  use  several  alphabets,  and  how  to  change 
from  one  to  another.  The  user  must  call  MXIALF  with  the  alphabet 
chosen  and  a  character  to  represent  it.  At  this  point,  any  time 
you  write  out  text,  this  character  causes  DISSPLA  to  switch  to 
the  alphabet  it  represents.  In  each  call  to  MESSAG,  these  special 
characters  are  not  printed  but  the  words  appear  in  different 
alphabets.  Each  string  must  end  in  the  BASALF  to  use  the  "$" 
string  counter  feature. 
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PROGRAM  EXAMPLE  15 
CALL  COMPRS 
CALL  SETDEV  (1.  1) 

CALL  BGNPL  (0) 

CALL  TRIPLX 

CALL  BASALF  ('STANDARD'} 

CALL  MX1ALF  ('STANDARD'.  '+') 

CALL  MX2ALF  (' L/CGREEK' ,  '#') 

CALL  MX3ALF  ('RUSSIAN'.  '/') 

CALL  MX4ALF  ('SCRIPT',  '&') 

CALL  MX5ALF  ('GREEK',  '$') 

CALL  MX6ALF  ('ITALIC'.  *!') 

CALL  SETEND  (';',  1) 

CALL  HEIGHT  (.5) 

CALL  PAGE  (l6. ,  8.5) 

CALL  PHYSOR  (0.,  0.) 

CALL  AREA2D  (l6.,  8.5) 

CALL  SETCLR  ('BLUE') 

CALL  MESSAG  ('THIS  IS  AN  1EXAMPLE+  USING  DIFFERENT;'.  100, 

•  5,  7.5) 

CALL  MESSAG  ('ALPHABETS.  SIT  IS  QUITE  EASY  !T0  USE  AND*; ' , 

100,  .5.  6.5) 

CALL  MESSAG  ('THE  ! ALPHABET  CAN  +BE  CHANGED  SAT  ANY  TIME*; ' , 

100.  .5.  5-5) 

CALL  MESSAG  ('NOTICE  1THAT  EACH  STRING  SENDS  IN  THE  +BASALF ; ' , 
100,  .5.  4.5) 

CALL  SETCLR  ('RED') 

CALL  MESSAG  ('SOME  MORE  EXAMPLES:;',  100,  2.,  3-5) 

CALL  MESSAG  ('RUSSIAN  -/RUSSIAN*  L/GREEK  -#ABCDEF+; ' , 

100,  .5,  2.5) 

CALL  MESSAG  ('GREEK  -  SABCDEFGH  +ITALIC  -1ITALIC+;', 

100,  .5,  1.5) 

CALL  SETCLR  ('BLUE') 

CALL  MESSAG  ('THE  ! STRINGS*  ARE  PLACED  IN  INCHES  FROM  ORIGIN;', 

100,  .1,  .5) 

999  CALL  ENDPL  (0) 

CALL  DONEPL 

CLOSE (1 .STATUS*' DELETE ') 

STOP 

END 


This  is  an  example  using  different 
Alphabets,  it  quite,  ooau  to  vise  and 
the  ALPHABET  can  be  changed  at  amu  time. 
Notice  that  each  string  m!  in,  the,  BASALF 
Some  more  examples: 

RUSSIAN  -PKXXHAH  L/GREEK  -a^pdcip 
GREEK  -  ABHAE$TX  ITALIC  -ITALIC 
The  strings  are  placed  in  inches  from  origin 
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Instruction  Alphabet 


This  next  example  shows  the  power  of  the  instruction 
alphabet.  The  user  can  choose  any  alphabet,  spacing,  height  or 
font  through  commands  similar  to  those  of  multiple  alphabets.  The 
commands  are  enclosed  in  parentheses  and  operate  only  on  the  line 
they  appear  in.  Below  are  the  commands  that  were  used  for  example 
16,  for  a  complete  list  refer  to  the  reference  manual  or  pate  B- 
24  of  the  Pocket  Guide. 

Pi  -  set  the  i-th  tab  at  the  current  position  horizontally.  (i<20) 
Gi  -  move  back  to  the  i-th  tab  previously  set  by  P.  (i<20) 

Er  -  move  up  (superscript)  from  base  line  a  distance  r. 
lr  -  move  down  (subscript)  from  base  line  a  distance  r. 

Yr  -  skew  character  by  factor  r*x/y. 

Ui  -  underline  once  from  i-th  P  tab  to  current  position. 

Di  -  double  underline  from  i-th  P  tab  to  current  position. 

Fi  -  Font  changed  to  corresponding  style: 


0  -  Default 

1  -  CARTOG 

2  -  SIMPLX 

3  -  SCMPLX 

4  -  COMPLX 

5  -  DUPLX 

6  -  TRIPLX 

7  -  GOTHIC 


8  -  FUTURA 

9  "  SERIF 

10  -  FASHON 

11  -  LOGOI 

12  -  SWI SSL 

13  -  SWISSM 

14  -  SWISSB 


NOTE:  If  a  shaded  font  is  used  it  must  be  called  prior  to  use. 
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$ 


Mi  -  Set.  character  set 

to  specified  alphabet: 

0  -  STANDARD 

7  -  L/CGREEK 

1  -  L/CSTO 

8  -  RUSSIAN 

2  -  ITALIC 

9  -  L/CRUSS 1  AN 

* 

3  -  L/C  ITALIC 

10  -  SPECIAL 

1*  -  SCRIPT 

11  -  MATHEMATIC 

5  -  L/CSCRIPT 

12  -  HEBREW 

6  -  GREEK 

13  -  ALFBET 

X  -  as  an  argument,  resets  instruction  to  default  value. 


Example  #16  writes  out  text  in  various  styles  using  the 
instruction  alphabet.  Text  can  be  slanted,  underlined  and  at 
different  heights.  Also,  the  mathematical  formula  is  written  out 
using  the  mathematical  character  set  (see  manual).  If  time  is 
taken  to  learn  the  instruction  alphabet,  the  user  can  benefit  by 
having  every  style  of  character  DISSPLA  can  produce  available  to 
him  throughout  his  plot. 
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PROGRAM  EXAMPLES 
CALL  TK41  (4105) 

CALL  SETDEV  (1,  1) 

CALL  BGNPL  (0) 

CALL  MIXALF  (' INSTRU') 

CALL  HEIGHT  (.2) 

CALL  PAGE  (11.,  8.5) 

CALL  PHYSOR  (0.,  0.) 

CALL  AREA2D  (11.0,  8.5) 

CALL  SETCLR  ('GREEN') 

CALL  MESSAG  ('  (H2F6M0C)  THIS  IS  AN(Y.5)  EXAMPLE  (YO.)  $' , 

.  100,2. ,7-0) 

CALL  MESSAG  ( '  (H2F7MO) OF  THE  (F6P1H3) INSTRUCTION  (H2U1)  (P2) 
ALPHABET  (02) $' ,  100,  1.,  6.) 

CALL  MESSAG  (' (H2F6M10C) E  (MO)-  MATHEMATICS  -  (M10) E$ ' , 

100,  2.,  5-) 

CALL  SETCLR  ('RED') 

CALL  MESSAG  ('  (M0F&H2P  1)  L I M  (G 1 L 1  .5H 1M1 1)  L  R  L  T  (M0)0(H3) 
(Mil)  (P2E . 2) S  (G2MOL2.5H1)  I  -  1  (G2LXE4 .5$3) N 
(EXM1 1H2S3) 2$ ' ,  100,  1.,  3-0) 

CALL  MESSAG  ('  (F6H3M1 1P3) I  (G3L2.5M6H1) L (L2-7H.5M0)  1 

(H1L2.5M10F6)  8  (M0)Y(L2.7H.7)  I  (L2.5M10H1)9 
(G3LXE3-5H1M6S7)L  (E3-3H.5MO)  2  (E3.5M10H1)8(M0)Y 
(E3 . 3H .5)  I  (E3-5H 1M10) 9$ ' .  100,  'ABUT',  'ABUT') 
CALL  MESSAG  ('  (F6EXMOH2)  F  (Ml  1)  8  (MO)  X  ,  Y  (L . 5H 1 F6)  I  (S2)  S  '  , 
100,  'ABUT',  'ABUT') 

CALL  MESSAG  ( '  (LXH2M1 1 F6) 9  (MO)  DX  (Ml 1H2)  3  (M 1 1 )  R  (M0L.5H1) I 
(LXH2) Y$ 1 ,  100,  'ABUT',  'ABUT') 

999  CALL  DONEPL 

CLOSE  (1 , STATUS- 'DELETE  1 ) 

STOP 

END 


NOTE:  in  the  call  to  MESSAG,  when  ABUT  is  given  as  the  position 
for  the  X  and  Y  value,  the  message  is  placed  where  the  previous 
message  left  off. 
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DTNSRDC  ISSUES  THREE  TYPES  OF  REPORTS: 


1  OTM8RDC  reports,  a  formal  aarlaa,  contain  information  of  permanent  technical 
value.  They  carry  a  consecutive  numerical  identification  regardless  of  their  classification 
or  the  originating  department. 

2  Departmental  reporta,  a  aemtformal  aeries,  contain  information  of  a  preliminary, 
temporary,  or  proprietary  nature  or  of  limited  interest  or  significance.  They  carry  a 
departmental  alphanumerical  identification. 

3.  Technical  memoranda,  an  informal  series,  contain  technical  documentation  of 
limited  use  and  interest.  They  are  primarily  working  papers  intended  for  internal  use. 
They  carry  an  identifying  number  which  indicates  their  type  and  the  numerical  code  of 
the  originating  department.  Any  distribution  outside  DTNSRDC  must  be  approved  by 
the  head  of  the  originating  department  on  a  case-by-case  basis. 


I 


OTNSflOC  5402/24  (R*v  4-M) 


/A-™> 

DTIC 


